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Abstract. We examine the effect induced on cosmological correlators by the simultaneous 
breaking of parity and of statistical isotropy. As an example of this, we compute the scalar- 
scalar, scalar-tensor, tensor-tensor and scalar-scalar-scalar cosmological correlators in pres¬ 
ence of the coupling C = \F 2 + \FF) between the inflaton cj) and a vector field with 

vacuum expectation value A. For a suitably chosen function /, the energy in the vector field 
PA does not decay during inflation. This results in nearly scale-invariant signatures of broken 
statistical isotropy and parity. Specifically, we find that the scalar-scalar correlator of primor¬ 
dial curvature perturbations includes a quadrupolar anisotropy, P^( k) = P(k)[l + g*(k- A) 2 ], 
and a (angle-averaged) scalar bispectrum that is a linear combination of the first 3 Legendre 
polynomials, B^(k\, fo, fci) = Y,l°lPl (ki ■ k 2 )P(fei)P(/c 2 ) + 2 perms, with co : ci : C 2 = 2 : 
—3 : 1 (ci / 0 is a consequence of parity violation, corresponding to the constant 7 / 0). 
The latter is one of the main results of this paper, which provides for the first time a clear 
example of an inflationary model where a non-negligible c\ contribution to the bispectrum 
is generated. The scalar-tensor and tensor-tensor correlators induce characteristic signatures 
in the Cosmic Microwave Background temperature anisotropies (T) and polarization (E/B 
modes); namely, non-diagonal contributions to {ae imi a^ 2m2 }, with \t\— £ 2 ! = 1 in TT, TE, EE 
and BB, and \(,\ — = 2 in TB and EB. The latest CMB bounds on the scalar observables 

(< 7 *, co, ci and C 2 ), translate into the upper limit pa/P< j> 1$ 10 -9 at 7 = 0. We find that the 
upper limit on the vector energy density becomes much more stringent as 7 grows. 
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1 Introduction 

Pseudo-scalar fields, which can naturally emerge in global symmetry breakings, have often 
been employed in models of cosmological inflation (e.g., see refs. [1-11]). A typical model of 
a pseudo-scalar coupled to a gauge field is 

£ = ~\ W) 1 2 - V (0) - \f 2 - ^FF , (1.1) 

where <fi is the pseudoscalar field and 1// expresses the strength of the axial coupling. Due 
to the motion of the inflaton 4>(t), this coupling enhances one of the two helicity states of 
the gauge field during inflation, inducing parity violation. This violation can be imprinted 
on the gravitational waves through the gravitational interactions of the gauge field. Most 
common observables of such chiral tensor perturbations are the cross correlations between 
temperature anisotropies T/E-mode polarization and B-rnode polarization (TB/EB) of the 
Cosmic Microwave Background (CMB) [3, 12-15]. Since the inflaton typically speeds up 
during inflation, the production of gauge quanta can potentially increase at smaller scales, 
giving rise to a gravity wave signal observable at terrestrial interferometers [16-18]. The 
chiral nature of this signal can be probed by combining measurements from multiple inter¬ 
ferometers [19, 20]. In addition, the f(cf>)FF interaction can also generate large primordial 
non-Gaussianity [13, 21-28], and primordial black holes [29]. Finally, this coupling has also 
been employed in models of inflationary magnetogenesis [30-38] . 

In the framework of the pseudoscalar inflation, impacts due to statistical anisotropy 
have been recently analyzed, by introducing a coherent vacuum expectation value (vev) of 
the gauge field [39-41]. In such a case, the primordial correlators are affected by the simulta¬ 
neous breaking of parity and statistical isotropy. This generates characteristic (ae imi a e n m 2 ) 
correlators of the CMB temperature anisotropies and polarization. Specifically, one finds 
non-vanishing TT, TE, EE and BB correlators for \l\ — £ 2 | = 1, and non-vanishing TB 
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and EB correlators for \t\ — = 2 [41]. Mathematically, all the CMB signals satisfying 

\l\ — h | = odd (even) in TT, TE, EE and BB (TE and EB) are parity-odd. These can be 
realized only when the primordial correlators include parity-violating information. On the 
other hand, the primordial correlators can generate nonzero CMB signals satisfying £\ ^ £ 2 , 
only when they break rotational invariance. Hence, the above signals become distinctive 
indicators of broken parity and rotational invariance in the primordial Universe. 1 

In ref. [41], we have presented a proof of the existence of these interesting signals, 
analyzing the simplest lagrangian (1.1). In this case, however, the vev of the vector field 
(which is the quantity breaking statistical isotropy) decays very rapidly during inflation, 
PA oc A 2 oc a -4 (where a is the scale factor), so that one needs to assume that a vector vev 
is present as an ad-hoc initial condition when the CMB modes left the horizon. The rapid 
decrease of p\ then implies that only the largest CMB modes can be affected by it, so that 
the induced non-diagonal correlators are highly red-tilted (they are oc AT 4 ). It was found in 
ref. [41] that only CMB multipoles with t < 10 can be affected at a detectable level. 

As anticipated in ref. [41], a more natural initial condition, and a more interesting signal, 
can be obtained if the kinetic term of the vector field is modified as in the Ratra mechanism 
[42], so to allow for a nearly constant p\- Specifically, if the kinetic term is — ^pF 2 , and if 
the functional form of /(</>) is chosen such that the background evolution satisfies = 

a~ 4 (t), then the vector has a constant electric vev (we are using standard electromagnetic 
notation for simplicity, although we do not need to assume that A M corresponds to the 
standard model photon). This time dependence can be achieved through a suitable relation 
between f(4>) and the inflation potential [43]. 

A model that provides a non-decaying pa and that breaks parity has been recently 
considered in ref. [36] as a model for primordial magnetogenesis [36]. It is characterized by 
the lagrangian 2 


z = -\ W) 2 - (~\f 2 + Iff) , ( 1 . 2 ) 

where 7 is a constant. 3 Analogously to the well studied 7 = 0 case [44-54], we show 
in this work that, for /(</>) oc a -2 , the interaction between 1 and A^ leaves nearly scale- 
invariant signatures on the primordial cosmological correlators. Differently to the 7 = 0 
case, however, due to the / 2 (0)FF interaction, one vector helicity state is produced with a 
greater abundance that the other one. As mentioned above, this violation of parity can affect 
the CMB correlators through the gravitational interactions of the vector field. 

In addition, a stronger effect arises from the direct coupling between the vector field and 
the inflaton. We mentioned that the interaction in eq. (1.2) is responsible for (i) maintaining 
a nearly constant energy in the vector field during inflation, and (ii) enhancing one helicity of 
the vector field with respect to the other one. These effects are due to the classical evolution 
of I(<f>) and therefore to the vev of the inflaton field. However the same interaction term 
(1.2) also couples the vector quanta to the inflaton perturbations. As we will see, in this 
model pa/P p <§C 1 is required. Therefore one can identify the perturbations of the inflaton 

1 These discussions rely on the assumption that there is no mechanism breaking parity and isotropy at late 
times. 

Tn the proposal of ref. [36], the field <f> does not need to be the inflaton field. 

3 The explicit breaking of parity can be avoided if the F 2 and the FF term are proportional to two different 
fields (as for instance in supergravity) that have vevs that evolve during inflation maintaining a constant ratio 

[36]. 
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as the adiabatic perturbation of the model, with negligible error. As we show below, the 
couplings 5<j)5A and 5(j)8A 2 encoded in eq. (1.2) modify both the power spectrum and the 
bispectrum of the primordial scalar perturbations. For the power spectrum of primordial 
curvature perturbations, one obtains the ACW [55] quadrupolar term 


Pc (k) 


P{k) l+ 5 *(k-A) s 


(1.3) 


with a nearly scale invariant g * parameter. For the (angle-averaged) bispectrum, one obtains 
the first three terms of an expansion in Legendre polynomials [52]: 

2 

B c (k 1 ,k 2 , k 3 ) = ^2 clPl {ki ■ k 2 )P(ki)P(k 2 ) + 2 perms., (1.4) 

L =o 


with co : ci : C2 = 2 : —3 : 1. 1 The three parameters { g *, co, C2} arise due to the breaking of 
statistical isotropy, and were already found in the model [51, 52]. The non-vanishing 

of ci instead requires both breaking of statistical isotropy and parity. This is the first time 
that, within these conditions, a non-negligible contribution ci to the curvature bispectrum is 
generated during inflation. 4 5 * 

As we show in section 4.1, for 7 = 0, the fact that such signatures have not been 
observed in the CMB data enforces the upper limit pa/ Pp ~ 1CT 9 . In our work we study 
the 7/O case. Our results are summarized in figure 1, where we show the upper limit on 
Pa/ P<p as a function of 7, starting from 7 > 1. Strictly speaking, our analytic computation is 
valid for |. However, the limit is continuous in 7, and figure 1 shows that extrapolating 
our lines in the 7 1 intervals provides results consistent with the limit at 7 = 0. To give 

a measure on how strongly the upper limit decreases with 7, we note that our computation 
provides the upper limit pa/ Pp ~ 4 x 10 -16 at 7 = 1, and pa/ Pp ~ 3 x 10~ 25 at 7 = 2. 

Although having such a small pa is a mathematical possibility, in section 4.2 we argue 
that a greater energy in the background vector field should be expected, simply from the 
addition of IR modes [51]. This essentially rules out the model (1.2), under the assumption 
that it produces a constant vector energy density, that 7 > 0(5), and that <fi is the inflaton. 
Our results do not immediately extend to the case in which is not the inflaton. However, 
we expect that also in that case the primordial perturbations will be significantly affected by 
the effects we have studied, mostly due to the linear coupling (which exists at least due to 
gravity) between 5<f> and the inflaton perturbations [25]. 

This paper is organized as follows. In the next section, we review an inflationary model 
based on the lagrangian (1.2) and we discuss how it can lead to a constant vev of the 
gauge field. In section 3, we compute the scalar-scalar, scalar-tensor, tensor-tensor and 
scalar-scalar-scalar correlators of primordial curvature perturbations, and find the distinct 
observable predictions of the lagrangian (1.2). In section 4, we estimate the observational 
bounds on the energy density of the gauge field vev from the latest Planck constraints on g *, 
co, ci and c 2 - This result is discussed in section 4 and in the concluding section. 

4 Notice that, due to a non-vanishing vev of the vector field, a statistical anisotropic bispectrum is actually 
generated and, after an angular average (see the discussion after eq. (4.2)), it takes the form (1.4). For studies 
of CMB bispectra that break statistical isotropy see refs. [56, 57]. 

5 See refs. [52, 58] for discussions on another possibility to generate non-negligible ci from large-scale helical 

magnetic fields at the radiation dominated era. 
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2 A model for breaking parity and statistical isotropy 


Let us consider the action [36] 6 


S= d A Xy/^g 


Ml 


1 


-Mr - -dMd»(j> - V((/>) + I\(t>) + -t-F^Ff 


1 


7 


fJ,V 


, ( 2 - 1 ) 


with the parameter 7 being a constant. Differently from ref. [36], we identify 4> with the 
inflaton field, and compute how the couplings in eq. (2.1) affect the primordial perturbations. 

For 7 7^ 0, the coupling of the inflaton Held to the vector explicitly breaks parity, as the 
product F 2 is a scalar, while FF is a pseudo-scalar quantity. This coupling also affects the 
inflaton perturbations, since expanding <f> = 4>o(t) + x), and denoting I(cpo) = Iq(t), 

one has 


I{4>) = h + ^rS(t> = Iq(t) + <5/( t , x ) . (2.2) 

In principle, 0(5(f> 2 ) terms in the expansion of /(</>) could be considered, and included in the 
computation of the bispectrum. As discussed in ref. [51] (see also ref. [59] for an explicit 
check), these higher terms give a subdominant contribution, and we can therefore disregard 
them. 

We assume that the function I evolves in time during inflation as 

7o(r) oc a n (r) , (2.3) 

and it then sets to a constant after inflation (when <p sets to a minimum). ' Without loss 
of generality, we can take this value to be 1. Strong coupling considerations put the n > 0 
regime into question [60]. The coupling of the vector to any field charged under the U(l) 
symmetry is oc j^. For n > 0, this implies a very large coupling during inflation. Even if no 
real quanta of the charged particles exist during inflation, loop of virtual charged particles 
are out of perturbative control, which puts in question any perturbative result obtained from 
the model [60]. 

2.1 Background evolution 

A non-vanishing vev of the vector field during inflation leads to anisotropic expansion. It 
is well known (see for instance ref. [51]) that, for 7 = 0, consistency with the CMB results 
requires that the energy density of the vector field is much smaller than that of <j>. We show 
below that the limit on the vector field energy becomes even stronger for 7 / 0. Therefore, 
we can neglect the departure of the background geometry from the FLRW metric, and we use 
ds 2 = a 2 (—dT 2 +dx 2 ) in our computations [51]. As in the standard case, inflation is supported 

by the inflaton potential, and the standard slow-roll condition applies e = ^ — 

6 We use the following notation: is the field strength of A^, while F' 11 ' = ^M=F a p is 

the dual tensor, with ? 7 ° 123 = 1. We also use M p = l/y/8nG, where G is Newton’s constant. In the following, 
dots (primes) denote derivatives with respect to physical (conformal) time, while H = -, where a is the scale 
factor. 

7 As we mentioned, the required time dependence can be obtained by suitably relating I with the inflaton 

potential. See ref. [43] for details. 
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HM V 


< 1, with negligible corrections. At zeroth order in slow roll, the dependence of 


the scale factor on conformal time is a ~ — ( Hr ) . 

We do not need to identify the vector field with the standard model photon, 
nonetheless use the “electromagnetic” convention with 


We 


e = _Mj) A '. b = ^ 


V x A 


(2.4) 


where the Coulomb gauge Aq = V • A = 0 has been assumed. In this notation, the energy 
density in the gauge field acquires the familiar form pa = E )> B ■ For convenience, we 
introduce the canonical field V = Io(r)A, and we expand it in a background value plus 
fluctuations as V = V^°^(r) + SV(t, x) (consistently with spatial homogeneity, we have 
imposed that the background value depends on time only; we note that with this choice the 
term proportional to 7 in eq. (2.1) does not contribute to the background evolution). The 
gauge field vev then satisfies 


- ^V (0) = 0 , 

Mi 


(2.5) 


leading to vanishing magnetic component, B (0) = 0, and to an electric component that 

I" 9 

depends on the parameter n. In particular, for n = —2 (or -^ = ^ 2 ), a time-independent vev 


arises as 


E (0) (r) = E vev 

This is the case that we study in this work. 


2.2 Gauge field fluctuations 

The fluctuations of the gauge field have two helicity states (A = ±1): 


5Vi(T, x) 


d 3 k 
(2-7t ) 3 / 2 


E '5if ) M4 A, (k)e ikx 

A=±l 


( 2 . 6 ) 


(2.7) 


where the quantum field SV is decomposed as 

Svi X \r) = a x (k)6V x (T,k) + a{(-k)SV^T,k) , (2.8) 

in terms of creation and annihilation operators that obey the algebra [a^(k), a^,(k')] = 
<5 aa'< 5® (k — k'). The quadratic action gives the evolution equation for the mode functions: 

SVj! + (V + 2A SV X = 0 . (2.9) 

The parity violating term results in a contribution that differs in sign for the two helicities. 
Following ref. [36], we define the coupling strength parameter 

£ = -" 7 , ( 2 - 10 ) 

8 The polarization vector e[ A, (k) satisfies feie[ A, (k) = 0, rf'^kie^ (k) = — Aie^fk), e[ A ^*(k) = e[ _A) (k) = 
e^ A) (—k) and e( A) (k)e( A , (k) = <5 A _v- 
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and for definiteness we assume £ > 0. As we are interested in n = — 2 , this means that the 
parameter 7 is positive. In the opposite case, one simply needs to interchange 5V+ «->• 5V- 
in the results below. Given that Iq oc a n , one finds, at zeroth order in slow roll, 

SV" + ^ k 2 + 2A*4 - SV X = 0 . (2.11) 

Ref. [36] provided the solution to this equation with the standard adiabatic vacuum 
initial condition for arbitrary n. They obtained a particularly simple expression for £ 1, 

and in the long wavelength regime. For n = —2, this expression reads 

^ 2 2 ^ 3/2 > |fcr|«^«l, ( 2 . 12 ) 

while the negative helicity mode SV- (r, k) is produced on a much smaller and negligible 
amount [36]. Using electromagnetic convention, we can write the following power spectra in 
the long wavelength regime 

(SXi(r, k )5Y j (r', k')> « 5X + (t , k)5Y + {r\ fe)e l (+) (k)ej +) (k , )5 (3) (k + k') , (2.13) 

where X. Y = E,B and 

p7r£ o jyt2 o 

S E + (T,k) = -^- 2 ^ j j J - 2 =-SB + (T,k). (2.14) 

This expression indicates that, on super horizon scales (—kr -C 1), the magnetic corre¬ 
lators are much smaller than the electric ones. For this reason, we may ignore the magnetic¬ 
mode contributions in the primordial curvature and tensor correlations computed in section 3. 

2.3 Curvature and tensor perturbations 

To compute metric perturbations, we decompose them into scalar and tensor perturbations 
as standard [61], and we work in the spatially flat gauge for the scalar sector, <5g?9 alar = 0. In 
this gauge, the spatial part of the metric fluctuation is given by the gravitational wave alone 
as 5gij = a 2 hij. Moreover, since the energy density in the vector field is much smaller than 
that in the inflaton, we can identify the scalar curvature perturbation £ = 
while the gauge field perturbations 5Ai are isocurvature modes that are produced as discussed 
in the previous subsection, and will affect £ and hij through their couplings to them. 9 

We have already discussed the quantization of the 5Ai modes. The curvature and tensor 
perturbations are decomposed, respectively, as 


c(T ' x) = /( 273/2^-. 

(2.15) 

Mux) ~ f (2?r)3/2 y, ^ ) w4 ) ( k ) elkx ’ 

A—i2 

(2.16) 


u In the spatially flat gauge, there are also metric perturbations Sgoo and Sgoi- These are non-dynamical 
modes that are integrated out, and induce gravitational couplings between the gauge quanta and As the 
term proportional to 7 is a topological term, no metric perturbations enter there, so these couplings originate 
from the I 2 F 2 term, and are given in ref. [62]. Such couplings are suppressed (technically, by an e factor [62]) 
with respect to the direct 5(f>5A 2 couplings present in eq. (2.1), and so we can simply disregard them, and set 
Sgoo = Sg 0i = 0 [41, 51, 62], 
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where we have used the polarization tensor given by a product of the polarization vector as 


. ( A ) ^ A ^ 

(k) = \[2eA (k)e^ 2J (k). The quantized fields are expressed with the operators satisfying 


[a A (k), at, (k 7 )] = <5 AA '<5^<5(k — k 7 ), respectively, as 


4(t) = a 0 (k)C(r,fc) + aJ(-k)C*(r,fe) , (2.17) 

X \t) = a A (k)/i(r, k) + a\(—k)h*(T, k) . (2.18) 


We treat the effect of the gauge perturbations perturbatively, and indicate the scalar 
and tensor correlators as 


( 


Cki Ck 2 


) 


(CkiCk^Ckb 


Cki Ck; 


1 SK 2 / o +^Ck 1 Ck 2/i 


Cki Ck 2 Ck 3 y + \^CkiCk 2 Ck 3/ ^ 

h^h^) +(hi Xl) tf Xa) 


ki ,l k 2 / 0 

(Aa)\ 


ki k 2 


Cki h 


k 2 J l 


(2.19) 


The suffix zero refers to the correlators at zeroth-order in the gauge fields, namely to the 
standard inflationary vacuum correlators. We disregard the zeroth-order scalar three point 
function, as it would correspond to the non-Gaussianity from standard single-field slow-roll 
inflation, which is unobservable [63-66]. 

In absence of gauge fields, we have the standard mode function solutions 


LA _ fc) _ iH{l + ikT) jkr 
U ’ ’ 2y/e 2^/eMpk 3 / 2 6 


if 5Ai = 0 , 


leading to the standard power spectra 


Cki Ck 2 


W 2 2) ) 0 


2 ^-V8^(k 1 +k 2 ) , V = 

^-eVS^i k 1 + k 2 )<5 Al , A2 . 


H 2 

8t T 2 eM 2 ’ 


( 2 . 20 ) 


( 2 . 21 ) 

( 2 . 22 ) 


The suffix 1 in eq. (2.19) denotes the first non-vanishing correction from the vector fields. 
This requires two interaction terms, that we compute at tree level in the in-in formalism (see, 
e.g. ref. [67]) in the next section. 


3 Correlators with broken parity and rotational symmetry 

In this section, we compute contributions from the gauge held to the correlators (2.19). As 
we mentioned, we treat the gauge held contributions perturbatively. The dominant coupling 
relevant for the scalar modes is obtained from the direct coupling I 2 (—jF 2 + ^FF). By 

expanding d(I 2 ((j))) — —2^^-(, this gives rise to a term oc £[7b°) + 5A] 2 (as we discussed 
above, higher order expansions of I 2 and couplings that originate from 6 goo and 5goi can be 
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neglected). From this expansion we have the two terms 10 


S C i = / drrx | — 


2 a 3 \ I' n 




S <2 = / <ir# X | -§-) |c 


) - 7 ( E^SBi + B^'SEi ) , (3.1) 


?(°); 


( 0 ); 




- (SEiSEi - SBiSBi) - -iSE^Bi 


(3.2) 


The I 2 F 2 term also leads to the dominant coupling between the gravity tensor mode 
and the gauge fields 


S hl = - f drd^aXj (E^SEj + B^SB^ , 


(3.3) 


where we disregard the h 5A 2 coupling as it does not contribute at tree level to two point 
correlators involving the graviton. 

In the constant vev case, i.e., 1$ oc a~ 2 = (— Hr ) 2 , these three terms give rise to the 
three interaction hamiltonians 


a 77>vev 

fffi(r) = - 1 1 ^ 


H 4 T 4 


d J p 5Ei ( t , p) £_ p (r) , 


Hp (r) = - 


2 f d 3 pd 3 p' 

H 4 r 4 J (2 tt) 3 / 2 


®i(hp)®i(hp')C-p-p'W , 


Hhi(j) = 


^vev 

H 4 t 4 


d 3 p6Ej(T,p)hij- p (r) . 


(3.4) 

(3.5) 

(3.6) 


In and we have dropped the contribution of 5Bi with respect to 5Ei, as the 
magnetic perturbation is much smaller than the electric one at super-horizon scales, see 
eq. (2.14) (the super-horizon regime dominates the time integrals of the in-in computation 
[51], as we discuss below). Since B^ = 0, and since FF oc E r B ,, this implies that the 
pseudo-scalar interaction does not contribute to the leading order expression of Hq\ and 
2 - For this reason, we recover the same dominant interaction hamiltonians as in ref. [51]. 
Nonetheless, the / 2 (</>) FF coupling strongly influences the phenomenological results, as it 
changes the gauge held mode functions, see eq. (2.14). 

In the following, we perform the explicit computations of the correlators using the in-in 
formalism. A useful intermediate result is the commutator between the zeroth-order fields. 
Using the zeroth-order mode functions (2.20), one finds 


Cki (rO.Ckafa) 


iH 2 
6 eM 2 

2 iH 2 

3 M 2 


[rf-r 3 ] <5 (3) (k 1 + k 2 ) , 

[ T 1 - T 2\ <5 (3) ( k l + ^Ai.Az • 


(3.7) 


3.1 Power spectra 


By means of the in-in formalism, the leading order correction due to the gauge held to the 
two point scalar correlation function is given by 


2 


n 4 „(r) 

77 .— 1 




n Ck„(T),#a( r i) 


n=l 




(3.8) 


1(, We note that the 2 £ tadpole is canceled once the exact equation of motion for £ is taken into account, 
and the solution (2.20) is properly modified. This effect is completely negligible, due to the fact that p^ w <C p 
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We are interested in this result at super horizon scales, — k n T <C 1. Modes of £ and of 5E 
appearing in this expression are the zeroth-order solutions, and so they are uncorrelated with 
each other. So, the expectation value in eq. (3.8) splits into two separate expectation values. 
The drt time integrals are dominated by modes in the long-wavelength regimes, since the 
small-wavelength modes are highly oscillatory, giving a highly suppressed contribution to the 
time integral [51]. On superhorizon scales, the electric field (contained in H^\ (r n )) becomes 
a classical (commuting) field [51]. Therefore, the expression (3.8) is evaluated to 


Ck„ 


\n=l 


4 E 


drl f Tl dr 2 , ,, 


H 4 j ./ rf J r; 


-f I d 6 Pi d 4 p 2 


Ej ev E] ev mirupJSEjfo, P2 )) 

Cki(r),C-pi(Ti)l) ([Ck 2 ( r ), C- P 2 ( T 2 ) ) + (ki o k 2 


(3.9) 


After computing this with eqs. (2.13) and (3.7), the time integrals are reduced to — f^ k -\ 77 , 


C T e drj 

Ti 

which is equivalent to the number of e-folds before the end of inflation at which the modes 
with k\ leave the horizon, given by N kl . With an identity: 


[ei +) (k 1 ) C W(-k 1 ) + eW(k 2 ) 6 $ +) (-k 2 )] ^ev^vev^ 3 ) (fcl + ^ 


e, (s) (ki)e5' s) (ki )EJ ev E] ev 6^ (k, + k 2 


,(-*) 


s=±l 


1 - ( ki • E vev ) 5^ (ki + k 2 ) , 


(3.10) 


we finally obtain 

(ncJ) 


\n= 1 


E 2 2n£ N 2 

^vev _«i_ 

2ne 2 M£ £ 3 k\ 


1 - (ki • E vev ) 6® (^ + k 2 ) , £ » 1 , (3.11) 


which indeed shows a non-vanishing quadrupolar asymmetry, and which is scale invariant up 
to the logarithmic dependence on k of N k . 

It is instructive to compare this result to that obtained for £ = 2y = 0 [51]. We find 

2 


Cki Ck 2 


Cki Ck 2 


47T £ 3 


5E 4 


5»i 


5=o 


5E 4 


5=o 


+ 5E_ 


5=o 


(3.12) 


In short the difference between our result (3.11) and the result obtained for 7 = £ = 0 is 
simply due to the difference between the wave functions of the sourcing gauge fields (the 
mode function in the present case is given in eq. (2.14); for £ = 0 one has instead 5E+ = 


5E_ = 


[51]). 


y/2k 3 / 2 

By an analogous computation, the leading gauge field contributions to the scalar-tensor 
and to the tensor-tensor correlators are found to be 




E 2 e 2n ^ 


neM£ £ 3 


tfX i)«< 3 >( k i + k 2 )S X ,,2 , (3.13) 
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The absence of the A = — 2 mode in these expressions indicates parity-violating correlators. 
Moreover, these correlators also break isotropy due to non-vanishing quadrupolar term. The 
gauge field contributions to the scalar-tensor and tensor-tensor correlator are suppressed, 
respectively, by a factor e and by a factor e 2 with respect to the contribution to the scalar- 
scalar correlator. An analogous behavior was also found in the F 2 + c f>FF model [41]. 


3.2 Scalar bispectrum 

Again using the in-in formalism, the leading order contribution of the gauge field to the scalar 
bispectrum is given by : 


nc k „(r)\ =' 3 

/ 1 \n=l 


\ (211) / 3 \ ( 121 ) / 3 \ ( 112 ) 

nc k » + n< k >) +(n4 n (r)) 


\n =1 

where we have defined 

/ 3 \ ( a ^ c ) 


\n=l 


\n= 1 


t-t . pt r T i r r 2 

II Ckr,.( T ) ) dn dr 2 / dr 3 


\n= 1 


n 4 n {r),H Ca (ri) 


n=l 


M C 6 (r 2 ) 




(3.15) 


(3.16) 


As for the power spectrum, the time integrals are dominated by modes in the super¬ 
horizon regime. Proceeding as in that computation, the first term in eq. (3.15) evaluates 
to 


( 211 ) 


n Ck„ 


= — i- 


\n= 1 


. 3 2 E 2 ev r dr 1 r Tl dT2 f T2 dr^ f d 3 pi^Pi 

H 12 J rf J J r| J (27r) 3 / 2 

E] ev El ev (SEifapJSEiiTup'JSEjfo, P2 )SE k (r 3 , p 3 )> 


d 3 p 2 / d 3 p 3 


Cki (t) i C—pi —( T l) 

+5 perms in k n 


Ck 2 (r),C-p 2 ('r 2 ) )( Ck 3 (t ), C—p 3 (js) 


(3.17) 


This expression evaluates to 
/ a \ ( 2n ) 
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\n=l 


Q TT'Z TTZ 
^-^vev 


27r 2 (2nf 2 e 3 M® 


<J (3) (ki + k 2 + k 3 ) e 47r « 
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US £*3 
^2 ^3 


el +) *(k 2 ) e ; +) (k 2 ) e l +) *(k 3 )ei +) (k 3 )E™Er 


1 1 
*1 ’ k 2 ’ 


dr\ 


r 3 - r 3 


7Max 

+2 perms in k n , 


fTl 

~^2 


dr 2 


r 3 - t 3 
T '2 


e> i. 


fTl 


dT 3 


- Tq 


(3.18) 


Namely, the internal time variables in the integrals cover the region t > T\ > r 2 , r 3 , 
plus lower bounds dictated by the requirement that all the modes in the interaction are in 
the super-horizon regime. Starting from the other two terms in eq. (3.15), and relabeling the 
internal times, we obtain an expression where the integrand is identical to the integrand of 
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eq. (3.18), but the internal times cover the two complementary regions r > r 2 > 73 , t\, and 
r > 73 > n, T 2 - This leads to 


n 


\n=l 


„ 3E 2 ev # 2 ^ (ki + k 2 + k 3 ) e 4 ^ 

, ~ 27r 2 (2vr) 3/2 e 3 M® & 2 £ 6 

4 +l *(k 2 ) £ < +) (k 2 ) e < +) *(k 3 )4 +) (k 3 )E™"Er 


J Max 


1 1 1 


r 3 -r 3 r r 3 — r| f T , r“ - r; 
dn-,— / dr 2 - z — / dr 3 


3 ^3 
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Ti 


ki’ k 2 ’ kg 

+2 perms in k n , £ 1 . 


i 

’ fc 2 


To 


1 

’ fc 3 


To 


(3.19) 


The time integrations are dominated by the earlier possible times, giving 


tt /- \ _ < 5 ( - 3 ' ) (ki + k 2 + k 3 ) 

J_J_Ck n ) — /n_\S/9 ^ 1 ,. l7. wvev/l^ 1 ’^ 2 ’ ^ 3 ) 0 


in=l / i 


(27t) 3 / 2 


ki,k 2 ,E v 


+2 perms in k n , £ i§> 1 . 


(3.20) 


We here have defined 


k 1 .k 2 ,E vev 

f(h,k 2 ,k 3 -, e 


^4 +) *(ki) e ' +) (k,)4 +), (k 2 )4 +) (k 2 )E-“Er 

_ 3£, 2 „7f 2 Min[iV tl ,JVt 2 ,A' fe ]iV tl iV fa 


27r 2 e 3 M® £ 6 


r.3 r.3 

A ,1 ft , 2 


e»i 


(3.21) 

(3.22) 


where we recall that is the number of e-fold before the end of inflation at which the mode 
with momentum ki left the horizon. Using the identity e- + ^*(k)ey + '(k) = ^ — kikj + i rf l i k kk 

the angle dependence in ^, vev can be simplified as 


•7?>i 

ki, k 2 , E v 


| 11 - (kr ■ E vev ) 2 - (k 2 • E vev ) 2 + (ki ■ E vev ) (k 2 • E vev ) (k, ■ k 2 


- ki • k 2 + (ki • E 


+ i 


k 2 • E vev 


E vev • ki - k 2 


E vev - f ki x 


k2 )]} 


(3.23) 


In the I 2 ((j))F 2 model (the 7 = £ = 0 case), the corresponding contribution to the 
bispectrum is [51] 


n ) 

\n=l / 1 


4=0 


' ( 3 )(k ^ 2 +k 3 ) Ck,^/<fa.fe'fe^-o) 

+2 perm in k n , 


(3.24) 


11 





















where 


d =0 

ki,k 2 :E vov 


f(ki,k 2 ,k 3 ] £ = 0) 


_si=±l 

1 - (kq • E vev ) 2 - (k 2 • E vev ) 2 

+ (ki • E vev ) (k 2 • E vev ) (k! • k 2 ) , (3.25) 


,S9=±1 


4" l (fe)4‘ 1, ’(fe) 


jjjvev jjjvev 

j k 


f (ki,k 2 ,k 3 ; £ > 1) x 




4 


(3.26) 


with 5E 


3 H 2 


^ ^ 2 k z /2 being the mode function of either helicity mode in the £ = 0 case 

(where the two helicities are produced in equal amount). As for the power spectrum, the 
difference between the £ 1 and the £ = 0 result is simply due to the difference of the 

sourcing gauge fields. In the present case, one helicity dominates the final result, leading to 
violation of parity. 


4 CMB constraints 

In this section, we compare the primordial correlators induced by the action (2.1) to the 
latest CMB data, and we obtain an upper bound on the energy density in the background 
gauge field, p v £ v . Finally, we compare this bound with the theoretically expected value for 
the vev. 

4.1 CMB observables and constraints 

The primordial correlators computed in section 3 act as the initial conditions of the CMB 
correlators and can therefore be measured by CMB observations. Formally, the signatures in 
the CMB scalar-scalar, scalar-tensor and tensor-tensor power spectra are essentially identical 
to those obtained in the F 2 + <pFF model [41], since the primordial correlators have the 
identical angular dependence. The main difference between the signatures on that model 
and the signatures that we are computing here is that in the present case the gauge field 
vev is constant, so that these signatures are nearly scale invariant. On the contrary, in the 
F 2 + cj)FF case, the vector vev (that needs to be assumed as an ad hoc initial condition 
present when the CMB modes were generated) is rapidly decreasing, leading to observational 
effects only at the largest scales. 

As proven in ref. [41], the scalar-scalar correlator (3.11) can create TT, TE and EE in 
\l\ — 1 2 \ = 0,2, while the scalar-tensor (3.13) or tensor-tensor (3.14) correlator can generate 
TT, TE, EE, BB, TB and EB in \h - l 2 \ = 0,1,2. Specifically, TT, TE, EE and BB (TB 
and EB) in \l\ — 1 2 \ = 1 {\i\ — 1 2 \ = 2) are distinct signatures of the anisotropic pseudoscalar 
inflation, since these appear only in the case that parity and rotational symmetries are broken 
at the same time. However, these signal-to-noise ratios are smaller than the scalar-scalar ones 
because of the slow-roll suppression of the tensor mode. Likewise, the tensor mode also gives 
special correlations in the bispectrum due to broken parity and rotational invariance; also 
such contributions are subdominant due to the smallness of the tensor mode. In this section, 
we therefore focus on the observables associated with the primordial scalar mode £. 
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For this analysis, we use the conventional g* parametrization for the power spectrum 
[55] and the cl parametrization for the bispectrum [52], namely, 


2 


n ^ 

n=l 


nc. 


n= 1 


<5® (ki + k 2 )P(fci) 

£ (3) (ki + k 2 + k 3 ) 
(27 t ) 3 / 2 


1 + 5* (ki • E VI 
^2 CLP{ki)P{k 2 )PLiXi ■ k 2 ) + 2 perms , 


(4.1) 

(4.2) 


where Pl{x ) is the Legendre polynomial. The latter bispectrum form is obtained after 
averaging the original anisotropic bispectrum (3.20) over all directions of E vev , in the spirit 
of isotropic CMB measurements. This is the quantity that is immediately associated to 
the angle-averaged reduced bispectrum computed from the data: the reduced bispectrum 
b (ki, k 2 , A&) is obtained by averaging the bispectrum over all possible orientation of triangles 
of sides of length k\, k 2 , k%. The theoretical prediction for the reduced bispectrum associated 
to eq. (3.20) is therefore equivalent to the theoretical prediction associated to the average of 
eq. (3.20) over all possible direction of E vev [52]. 11 

Keeping into account that the CMB data force |<7*| <C 1 [68, 69], the sum of eqs. (2.21) 
and (3.11) yields P(k) ~ and 


£»1 : 


9* - ~ 


1 9 /V 2 p 2?r £ n ve ' 

^ 'CMB e _ Pe_ 

ire £ 3 p ^ 


(4.3) 


where JVcmb is the number of e-folds before the end of inflation at which the CMB modes 
leave the horizon, p$ ~ V(4>) — 3M 2 H 2 is the energy density of the inflaton, and 
is the energy density of the gauge field. 

Let us now compute the average of the bispectrum (3.20) over all directions of E vev . 
Using 12 


d 2 F, v 


4-7 T 


CLe... = j) Po(kl ' + ^(k, • k 2 ) , (4.4) 


and setting — Nk 2 — Nk 3 — Nqmb, we obtain 


£ > 1 : c 0 = - 


4ATcmb e 2 ^ 

37t £ 3 


-g* 


ci = -- 


3c 0 


co 

C2 = ¥- 


(4.5) 


We compare eqs. (4.3) and (4.5) with the results obtained for the 7 = £ = 0 case [52] 

48A r pv TR p v p v 16 co . 

£ = 0 : 3* ----, c 0 = — -N C mb 9* , ci = 0 , c 2 = — . (4.6) 

e p</, 3 2 


The Planck collaboration reported the 95% CL limits [28, 69] 13 


- 0.0225 <g*< 0.0363 , -10.7 < c 0 < 16.7 , -89 < a < 324 , -57 < c 2 < 47 . (4.7) 

11 In the study of CMB anisotropic bispectra [56, 57] this would correspond to single-out a monopole term 
in a spherical harmonic expansion of the anisotropic bispectrum [56] that does contribute to the isotropic 
(angle-averaged) bispectrum. 

12 This identity is easily derived using the spherical-harmonics representations of a unit vector and polar¬ 
ization vector, see [70]. 

lli These cl bounds correspond to the temperature only limits [28]: = 2.5 ±5.7, /nl" 1 = —49 ±43 and 

/nl 2 = 0.5 ± 2.7 (68% CL), where c 0 = (6/5)/^ al , ci = -(12/5) and c 2 = -(48/5)/Nf; 2 hold. 
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Figure 1. Upper bounds on PjSfV/ty reported in eqs. (4.8) and (4.9), with IVcmb = 60 and e = 0.01. 
As all these bounds must be satisfied, the final bound on p™ v /P>£ is the most stringent of the bounds 
shown. Strictly speaking, the solid and dashed lines are based on an approximation that holds for 
|y| | (that is, |£| 1). We extended them up to |y| = 1, as a naive extrapolation of the lines 
shown in the figure is roughly consistent with the values for |y| = 0, that are exact. 


From these constraints we get rough estimates of the upper bounds on p™ v /using eqs. (4.3), 
(4.5) and (4.6): 


and 


7 = 0 


PT 

P<j> 


< 


1,3 x 10 9 (oat) (wsw) 
3.0 x 10 9 (gjy 

(obi) 


1.7 x 10 “ 8 


(from p*) , 
(from c 0 ) , 
(from c 2 ) , 


7»2 


pT 

p<t> 


< 


3.1 x 10 ~ 5 

1.1 x 10~ 4 


3 \ 2 

7 3 
e 4 ^ 

3 \ 2 
7 3 


(from p*) 


1.3xl0-(*)(^)(^ a 

(oat) (]vSk ) 3 ( fromc o) 

(oAT) (nsL ) 3 (from cr) 


1.7 x 10 4 (^y) (oat) (jvSi) (from c 2 ) 


(4.8) 


(4.9) 


We show these bounds in figure 1, for the special choices of IVcmb = 60 and e = 0.01. 
While at 7 = 0 the strongest bound on Pg v /p</, is obtained from the power spectrum, at 
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7> ^ the strongest bound is due to the bispectrum coefficient Co- This is due the fact that 
the gauge production increases with increasing 7, and the bispectrum is affected more than 
the power spectrum by this growth (at the technical level, the source for the power spectrum 
is oc 5E 2 , while the one for the bispectrum is oc 5E 4 . This generates the e 27r ^/£ 3 factor in 
co/ff* in eq. (4.5)). 

The result we have obtained assumed £ = 2y > 0, but it can be readily applied also 
to the case in which £ = 27 < 0. As clear from eq. (2.11), changing sign of £ = 2q simply 
amounts in interchanging the role of the two helicities of the gauge field. In the computation 
of the inflation correlators, one simply needs to replace £ —> |£| and interchange 0 e'- \ 
One can then verify that the limits (4.9) can be extended to the 7 <C — ^ region by simply 
replacing 7 with I7I. 

4.2 Comparison with the expected vector vev 

In this subsection we compare the limits on /p<p obtained in the previous subsection with 

what is naturally expected in the model (2.1). 

Beside producing scale invariant vector perturbations, the I oc a~ 2 case also sustains a 
constant vector vev. This is not accidental. Due to scale invariance, super horizon modes 
have the same power at any given time t. The modes also had the same power when they 
left the horizon. 11 Equal power at horizon crossing (that happens at different time for the 
different modes) and at the later common time t is obtained because the modes have frozen 
amplitude while outside the horizon. In the super-horizon regime a mode solves the same 
equation of motion as a constant vev. So, if the super-horizon modes are frozen outside 
the horizon, also a completely homogeneous vev is. This explains why the / oc a -2 choice 
produces both scale invariant electric perturbations, and a constant electric vev. 

In presence of the vector vev, the interaction between the gauge and the inflaton per¬ 
turbations leads to anisotropic correlators of the primordial perturbations. The fact that this 

anisotropy is not observed in the data results in a limit on the vev of the vector field, and on 

E 2 

the associated energy density p^ v = —3p-. The current limit is summarized in figure 1. 

Naively, one may simply assume that the gauge field has no homogeneous vev, in which 
case the limit computed here does not apply. While this is a mathematical possibility, such 
an assumption is very unlikely in a model that is constructed to give scale invariant gauge 
perturbations. Unless one makes the very ad-hoc assumption that the duration of inflation is 
limited to the one necessary to produce the CMB modes that we observe at the largest scales 
(that is, Wot = Wimb — 60), also perturbations of wavelength larger than our horizon were 
produced during inflation with a scale invariant spectrum. Each of these infrared (IR) modes 
is observed as a classical and homogeneous vector field on our sky [51]. We can actually 
observe only the sum of such modes, and this precisely constitutes the homogeneous vector 
vev E vev that we have studied in this work. 

If we could observe many disconnected patches of the universe that have a size compa¬ 
rable with our Hubble horizon, we would observe many realizations of the IR gauge modes 
generated during the first Wot — W?mb e-folds of inflation, and we would observe a differ¬ 
ent sum E vev of these modes in each Hubble patch. Under the hypothesis of Gaussianity 
of the gauge perturbations (even a 0(1) departure from Gaussianity would not change the 

14 This is just a statement of time translation invariance during inflation; this invariance is broken by the 
slow roll parameters, as the inflaton typically speeds up during inflation. This leads to departure from scale 
invariance - typical in slow roll inflation - that we are disregarding in the present discussion. 
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order of magnitude estimate performed here), the values of these sum would have a Gaussian 
distribution of zero mean and variance 


l7l » j : 


<E 2 (x)) = 


2 7r 2 


MR 


dkk 2 \8E + (k)\ 2 = 


9 H 4 e 47r l' y l 

647r 3 |7| 3 


f dk 

Itr k 


4 „47r|7| 


9 H 


647T 3 I'y | 3 


(Wot - Nqmb) 


(4.10) 


We observe only one Hubble patch, and hence only one of such realizations. The variance 
we have just computed gives the typical value of E 2 ev obtain in any realization. We therefore 
obtain 
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l7l » j ^ 


Pe 


P(j) 
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H \ 2 e 4?r l 7 l 


128x3 \WJ ITT Wo. - JVcmb) • 


(4.11) 


The corresponding values of g* and cl are 
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p/ ItI £ 
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10 


1 
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(4.12) 


Already requiring that the energy in the vector field is subdominant poses a significant 
limit. For instance, for |y| = 5.5, and for only 10 e-folds of inflation more than the last Womb 
ones (this is a conservative assumption, as typical models of inflation give a much longer 
duration), one finds a subdominant vector field only if H £ 10“ 13 M p (= 3 x 10' 5 GeV). The 
limit obtained in this work is much stronger. Combining the results (4.9) and (4.11), we find 
the bound 


1 

^ 2 


i£ ~ ° “ 


|-y| 9 / 2 / e / 60 \ 3 10 

e 6?r M V 0.01 V^cmb/ Wot — Womb 


(4.13) 


which, for |q| = 5.5, gives the bound jj- <, 10~ 43 . In principle, one can assume values of jj- 
extremely small during inflation, although this is clearly a challenge for model building. The 
absolute minimum is obtained by imposing an instantaneous thermalization of the inflaton 
decay products, with the minimum reheating temperature T r h. m i n — 4MeV required for BBN 


[71]. This gives jj- > 10 42 . We thus see that values |q| > 5.5 are ruled out. 
For comparison, in the 7 = 0 case one finds [51] 


7 = 0 


VK\ 

Pe 


P<j> 


typical 

H 


3 i7 2 

= Mj " JVomb) 


=> F i 1»10 


4-5 


I e / 60 y 10 

0.01 V WoMB/ Wot - WOMB 


(4.14) 


so that a much greater value of H during inflation can be tolerated in this case. 

The bounds (4.13) and (4.14) are shown in figure 2 for the specific choice of the infla¬ 
tionary parameters e = 0.01, Womb = 60, Wot = 70. Strictly speaking, the result (4.13) is 
valid only at |7| \ (that is, |£| 3> 1). However, the figure shows that a naive extrapolation 

of the line from the I7I S> 1 regime to the |7| < 1 one is qualitatively consistent with the 
result at 7 = 0. 
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Figure 2. Upper bounds on H reported in eqs. (4.13) and (4.14), with IVcmb = 60, e = 0.01 and 
IVtot — -/Vcmb = 10- One can see that values |q| > 5.5 are ruled out as they would require a reheating 
temperature smaller than that required for successful BBN. 

5 Discussion 

In ref. [41] we computed the non-diagonal correlators between multipoles of the CMB tem¬ 
perature anisotropy and polarization induced by the simultaneous breaking of parity and 
rotational invariance during inflation. A very simple way to realize such a breaking is to 
couple the inflaton cj) with some gauge field with non-vanishing vev through a pseudo-scalar 
interaction <f>FF. A vector field with a standard kinetic term is rapidly diluted away by the 
expansion of the universe. To cope with this, in the present work we considered the recent 
modification by Caprini and Sorbo [36] of the Ratra mechanism [42] for the generation of 
a primordial magnetic field during inflation. This model is characterized by the interac¬ 
tion / 2 (</>)(— jF 2 + | FF ), which explicitly breaks parity. For the suitable time dependence 
/((/>(t)) oc a 2 (t), where a is the scale factor, the model produces scale invariant magnetic 
perturbations, which is at the core of the mechanisms of ref. [42], For / oc a -2 the model 
produces scale invariant electric perturbations. This was first used in the 7 = 0 case to 
sustain anisotropic inflation [72], exploiting an “electromagnetic duality” in the model for 
/ j [73]. 

We restricted our computations to the / oc a -2 case; however, due to the duality, 
we expect that our results can be readily extended to I oc a 2 . We also assumed that the 
field cj) is the inflaton (this is not required in the magnetogenesis applications [36, 42]). In 
presence of the vector vev, the interaction between the gauge and the inflaton perturbations 


17 - 








leads to anisotropic correlators of the primordial perturbations. In particular, among these 
correlators, we have shown for the first time that it is possible to generate a primordial non- 
Gaussian signature during inflation proportional to c\ in the parametrization (1.4), which is 
due to broken parity and rotational invariance. 

The fact that the anisotropic signatures predicted by this model are not observed in the 

data results in a limit on the vev of the vector field, and on the associated energy density 

E 2 

Pg v = — S 82 -. The limit is already strong at 7 = 0, where standard values of the slow roll 
parameter e = 0.01 and of the number of observable e-folds of inflation Wcmb = 60 lead 
to Pe v /P< t> ~ 10~ 9 . Such a strong limit is due to the fact that the vector modes act as 
isocurvature modes that continue to source the inflaton adiabatic perturbations in the super¬ 
horizon regime, inducing a (^cmb) enhancement of the anisotropic two (three) point 

correlators, and to the fact that the direct inflaton-vector field coupling is stronger than the 
gravitational one [51]. The limit becomes much more stringent with growing | 7 |, see figure 
1 , as the amount of gauge quanta produced by the moving inflaton grows exponentially with 
|7l at lar S e I'Tl- 

The relevant question to ask is whether such stringent limits are compatible with the 
vev that should naturally be expected in the model. We have considered a mechanism 
that, by construction, gives scale invariant perturbations for the gauge field. In typical 
inflationary models, the primordial perturbations of size comparable to our horizon were 
produced IVcmb — 60 e-folds before the end of inflation. This does not mean that the scale 
invariance induced in the mechanism should stop at such wavelengths, but it is natural to 
assume that also modes produced earlier were scale invariant. 15 We denote such modes as 
IR modes. Each of these modes is observed in our sky as classical and homogeneous vector 
field on our sky [51]. We can actually observe only the sum of such modes, and this precisely 
constitutes the homogeneous vector vev E vev from the IR modes (that eventually sum up to 
a vev of the vector field from the classical equations of motion) that we have studied in this 
work. 

Let us denote with Wot the total number of e-folds of inflation, and assume that E vev = 0 
at the onset of inflation (if a sizable vev is already present at the start of inflation, we 
simply expect an even stronger effect, and even stronger bounds). The IR modes that 
are seen as homogeneous from the point of view of the CMB modes are those produced 
in the first Wot — W;mb number of e-folds of inflation. The value of the same on our 
Hubble patch is obtained as a stochastic addition, obeying the typical random walk relation 
(E^ ev ) | ex p ec t e d oc Wot — W:mb (this is completely analogous [51] to how condensates of light 
scalar field develop during inflation). This was studied in ref. [51] for the 7 = 0 case. In 
the present work, we extended this study to non-vanishing and large 7 , where, due to the 
significant increase of the amplitude of the produced gauge quanta, the limits imposed by this 
effect become much more stringent, see figure 2. For instance, by simply allowing inflation 
to last 10 more e-folds than the amount necessary to produce the CMB modes (this is a very 
conservative assumption), we could rule out I 7 I > 5.5. 

Our results hold under the assumption that the function I((j>) produces a constant 
electric field. This can be achieved for I(4>) oc a n , with n = —2. Using an “electromagnetic 
duality” of the model, this result can be immediately translated to the I oc a 2 case, although 
for n > 0 the mechanism is plagued by a strong coupling problem during inflation [60] (see 

15 One can make ad-hoc assumptions to prevent this: for instance one can postulate that the duration of 
inflation is only limited to IVcmb e-folds, or that the scaling I oc a ~ 2 did not hold at N > IVcmb- We do not 
make such ad-hoc assumptions in discussing what should be naturally expected from this mechanism. 
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the discussion after eq. (2.3)). On the other hand, the n/- 2 case remains to be studied. 
Moreover, it would be interesting to study how the limits obtained here weaken in the case 
in which the field cf> is not the inflaton. We hope to come back to these issues in a separate 
work. 
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